# Imports
import math as m
import numpy as np
import matplotlib.pyplot as plt
import holoviews as hv
import pylab as pl
from pylab import exp,cos,sin,pi,tan, pi
import pandas as pd
import seaborn as sb
import holoviews as hv
from IPython.display import SVG
import io
from PIL import Image
from random import random
import elastica as el
import elastica_neurons as en
from dynamics import *
%matplotlib inline
hv.notebook_extension()
m = 11
n = 11
tau = 6
# Number of orientation selective neurons
nosn = 9
timesteps = 200
# for vonMises function
k = 0.25
A = 3
np.random.seed(67)
ac_orient4 = np.pi*np.random.rand(m,n)
ac_orient4[2:8,6] = np.random.rand(6)/5
oneColor()
visualField(ac_orient4)
results, rs, direction, magnitude = runExperiment('el',m,n,nosn,ac_orient4,timesteps,tau,vis=True,el_factor=0.002,E0=4)
results
t = np.arange(0,timesteps,1)
plotOneField((0,0),t,nosn,rs)
ac_orient2 = np.pi/2*np.ones((m,n))
ac_orient2[5,5] = 0
visualField(ac_orient2)
timesteps = 400
results, rs, direction, magnitude = runExperiment('el',m,n,nosn,ac_orient2,timesteps,tau,vis=True,el_factor=0.0001)
results
t = np.arange(0,timesteps,1)
plotOneField((0,0),t,nosn,rs)
ac_orient3 = np.pi/2*np.ones((m,n))
for i in range(m):
ac_orient3[i,i] = 3*np.pi/4
visualField(ac_orient3)
results, rs, direction, magnitude = runExperiment('el',m,n,nosn,ac_orient3,timesteps,tau,vis=True,el_factor=0.0001)
t = np.arange(0,timesteps,1)
plotOneField((0,0),t,nosn,rs)
results